Hyperbolic Interfaces 



Luca Giomi 

School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA 

(Dated: September 26, 2012) 

Fluid interfaces, such as soap films, liquid droplets or lipid membranes, are known to give rise 
to several special geometries, whose complexity and beauty continue to fascinate us, as observers 
of the natural world, and challenge us as scientists. Here I show that a special class of surfaces of 
constant negative Gaussian curvature can be obtained in fluid interfaces equipped with an orien- 
tational ordered phase. These arise in various soft and biological materials, such as nematic liquid 
crystals, cytoskeletal assemblies, or hexatic colloidal suspensions. The purely hyperbolic morphol- 
ogy originates from the competition between surface tension, that reduces the area of the interface 
at the expense of increasing its Gaussian curvature, and the orientational elasticity of the ordered 
phase, that in turn suffers for the distortion induced by the underlying curvature. 



Understanding the origin of shape in structures at var- 
ious length scales is a central goal across many areas of 
science and engineering. Often, the shape is dictated by 
the balance of competing forces simultaneously acting on 
a system, or by the interplay between spontaneous geom- 
etry and the residual stresses accumulated during fabri- 
cation and growth [1 . A standard problem consists of 
finding the optimal shape of a system given its mechan- 
ical energy and the constraints the system is subject to. 
Conversely, given a target shape with specific geometri- 
cal properties, one could ask what is the simplest physical 
system in which this will occur as a result of mechanical 
equilibrium and kinetics. Such an inverse problem repre- 
sents a crucial step in several aspects of modern technol- 
ogy, ranging from robotics [2 to medicine [3 . In the case 
of soft elastic solids, such as thin sheets of polymer gels, 
controlling the amount of local growth or swelling was 
demonstrated to be a powerful technique to produce var- 
ious three-dimensional shapes of both positive and nega- 
tive Gaussian curvature [4 . Fluid interfaces also bear a 
great potential for designing surfaces with special or tun- 
able geometric properties. Soap films, droplets, vesicles, 
micelles and membranes are examples of this potential as 
well as of the intrinsic beauty of these fragile objects. 

When endowed with in-plane orientational order, the 
richness of shapes that can be attained by fluid inter- 
faces becomes enormous. In a recent article Xing et al. 
[5] have shown that nematic and smectic vesicles formed 
from block copolymers with liquid-crystalline side chains, 
self-assemble in a great variety of morphologies due to the 
competition between in-plane liquid-crystalline order and 
bending elasticity. In a spectacular experimental work, 
Gibaud et al [6] have reported the formation of several 
complex structures in colloidal membranes, arising from 
the interplay between chiral order and interfacial tension. 
The morphogenetic force in these systems originates from 
the fact that orientational order is frustrated by the Gaus- 
sian curvature of the underlying surface [THS] • When par- 
allel transported on a surface of non-zero Gaussian curva- 
ture, a tangent vector rotates with respect to its original 
orientation and the amount of rotation is proportional to 



the Gaussian curvature experienced along the path. For 
this reason, local orientational order, whether nematic, 
smectic or chiral, is inevitably hindered in presence of 
Gaussian curvature. 

In this Letter I show that a special class of surfaces of 
constant negative Gaussian curvature, can be obtained 
in fluid interfaces equipped with an orientational ordered 
phase. This remarkable phenomenon originates from the 
competition between surface tension, that reduces the 
area of the interface at the expense of increasing its Gaus- 
sian curvature, and the orientational elasticity of the or- 
dered phase, that in turn suffers for the distortion in- 
duced by the underlying Gaussian curvature. Hints of 
this behavior were observed by Frank and Kardar ^10] . 
who reported an example of defect-induced buckling of 
a nematic membrane into the pseudosphere. The latter 
can be considered a particular manifestation of a more 
general phenomenon that, as we will see in the follow- 
ing, does not require defects nor any specific rotational 
symmetry of the ordered phase. 

Let us consider a two-dimensional interface M of sur- 
face tension a, equipped with an orient ationally ordered 
phase. This can arise in various soft materials, such as 
nematic liquid crystals, tilted molecular chains, or hex- 
atic suspensions. Orientational order on a surface can 
be described via a unit vector n in the tangent plane of 
the surface. The total energy of the system can then be 
expressed in the form: 

E = j dA{(j^K\Vn\^) , (1) 

where V indicates the covariant derivative in the metric 
of M, and n is an orientational stiffness constant with di- 
mensions of energy. The orientational contribution in Eq. 
([T]) describes, the one elastic constant approximation, a 
generic p— atic phase, where local orientations are defined 
modulo 27r/p [7 . p = 1 represents polar order, p = 2 
corresponds to a nematic phase and p = 6 characterizes 
hexatic interfaces. The energy has a characteristic length 
scale, in addition to that set by the boundary, given by 
£ = y^Kja. Thus we might expect the bulk Gaussian cur- 
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vature of the interface to be of order K = — with the 
negative sign due to the natural hyperbolicity of fluid in- 
terfaces. Remarkably, this simple estimate is exact within 
a range of material and geometrical parameters, resulting 
in an interface of constant negative Gaussian curvature, 
whose intrinsic geometry is that of the hyperbolic plane. 
The remainder of this paper will be devoted to prove this 
result and explore its consequences in the special case of 
tubular interfaces. 

As a starting point, we express the vector field n in a 
local orthonormal frame of tangent vectors Cq, (a = 1, 2), 
so that: n = cos^ei + sin ^62- Using standard manipu- 
lations 1^ , the gradient terms in Eq. ([T]) can be then 
expressed as Vin = {diO — Ai) where Ai = ei -9^62 
is the spin connection that accounts for the rotation of 
the local frame as we move along the surface and 
n_\_ = cos ^62 — sinOei. The Gaussian curvature of 
the surface is given by the curl of the spin connection: 
K = e^^ViAj^ with e*-^ the antisymmetric Levi-Civita 
tensor. As a consequence, the orientational contribution 
to the energy ([T]) cannot vanish everywhere because A 
cannot equate the curl-free vector field VO on a surface 
having non-zero Gaussian curvature [8 . Orientational 
order is therefore geometrically frustrated. 

Minimization of Eq. ([T]) with respect to the orientation 
field 6>, leads to the well known equation [12]: 



form: 



V • (Vi9 - A) = , 



subject to the constraint: 



e'^Vi{djO-Aj)=r]-K. 



(2) 



(3) 



where r] is the topological charge density, should topolog- 
ical defects be present in the configuration of n [7 . The 
orientation at the boundary dM can either be prescribed, 
in which case 6 = 0{s) at dM (with s the arc- length of 
the boundary curve), or free. In the latter case energy 
minimization further demands ly • (V6> — A) = at the 
boundary. Here ly is the outward directed tangent vector 
normal to the boundary curve dM. A classic strategy to 
solve Eq. ([2| under the constraint ([3|, is to express the 
vector field — A through the antisymmetrized deriva- 
tives of a scalar field (p: diO — Ai = e^dk^P- With this 
choice, Eq. ([2| is automatically satisfied, while Eq. ([3| 
implies: 



Ag(p = r] - K 



(4) 



where = V^V* is the Laplace-Beltrami operator in 
the metric of M. As it was observed by Vitelli and Nel- 
son the scalar field (p plays the role of a geomet- 
ric potential expressing the amount of non- uniformity in 
the vector field n caused by the presence of topological 
charge, whether localized in the defects or distributed in 
the form of Gaussian curvature. In terms of the en- 
ergy ([1]) of the interface can be rewritten in the simple 
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= J dA{(j- 



(5) 



In the case of free orientation at the boundary, t • Vip = 
—U'{\/0—A) = 0, hence ip is constant along the boundary 
curve and can, in particular, be set to zero as the energy 
([5| only depends on the gradient of (p. Analogously, if n 
is everywhere tangent to 9M, ly • \/(p = t • (V^ — A) = 
f<ig^ where f<ig is the geodesic curvature of the boundary 
[13j . With the help of the Gauss-Bonnet theorem [13l 
[T4] . however, one can prove that this latter condition is 
automatically satisfied by any solution of Eq. ([5|, thus 
we can again take = at dM. 

Next, we calculate the energy variation corresponding 
to a displacement of the interface along the normal di- 
rection: r ^ r + eA/", with R the position vector and e 
a small displacement along the normal vector N. The 
variation of the area element is a classic result of differ- 
ential geometry: S{dA) = —2HedA with H the surface 
mean curvature [14]. In order to calculate the variation 
of the orientation energy, we make the assumption that 
the system is defect free, so that r] = and the second 
term in ^ can be expressed as: 

W = J dA\\/ip\^ = - J dAdA' G{r,r')K(r)K(r') , 

where G(r, r') is the Laplacian Green function, with 
G(r, •) = G(-,r) = for r G dM. Next, we intro- 
duce a system of conformal coordinates {u^^u^)^ so that 
the metric is locally Euclidean up to a space dependent 
scaling factor: ds'^ = w[{du^Y ^ {du^Y] |7]. In conformal 
coordinates, the area element is given by dA = w dv}dv? 
and the Gaussian curvature can be expressed in the sim- 
ple form: 
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A \ogw 
2w 



(6) 



where A = w/S.g = + is the Euclidean Laplacian. 
In the new coordinates the functional W becomes: 

W = -^J d^u A\og w{u) J d^u'G{u,u')A\ogw{u'), 

where d'^u = du^dv?. With these simplifications, all 
the geometrical structure of the interface is embodied 
in the conformal weight w and the expression of W is 
now suitable to calculate the normal variation. Consis- 
tently with the area variation given above, one has that 
5w = —2wHe and d{logw) = —2He. Thus, after some 
manipulations, we find: 

SW = J d^uA{He) J d^u'G{u,u')A\ogw{u'). 

Now, the integral over u' is a function of u that vanishes 
when u G dM. Thus, using Green's identities and taking 



3 




2h 



FIG. 1: A portion of catenoid (orange) superimposed by a 
portion of the pseudospherical hyperboloid described in text 
(green) . Both surface are bounded by a pair of rims of radius 
R held at a distance 2h from each other. 



into account that e = at the boundary, we find: 

SW = J d^uHeA\ogw{u) = j dAKHe . 

Combining this with the area variation we finahy obtain 
an expression for the normal variation of the total energy: 



6E = j dA {cr ^ nK){-2He) 



Evidently, the variational equation 5E = has two so- 
lutions: H = and, as anticipated, K = —ajn. This 
implies that, in order for the interface to have minimal 
energy, it can either form a surface of zero mean cur- 
vature, such as a classic soap film, or a surface of con- 
stant negative Gaussian curvature. In presence of topo- 
logical defects, the Gaussian curvature in Eq. ([7|) is re- 
placed K — T] and energy is again minimized by setting 
K = —ajn outside of the defect-core [10^. As we ob- 
served in the introduction, this behavior originates from 
the competition between surface tension, that acts to re- 
duce the area of the interface with a consequent increase 
in curvature, and the orient ational elasticity of the or- 
der phase, that in turn is subject to a stronger distortion 
as the Gaussian curvature is increased. For this reason, 
we might also expect that the stability of the purely hy- 
perbolic shape would be lost for small values of when 
the distortion induced by the curvature of the interface 
is energetically inexpensive. In the remainder of the pa- 
per, we will consider the example of a tubular interface 
to clarify these assertions. 

Let us consider an interface bounded by two circular 
rims of radius R held at distance 2h from each other (Fig. 
[T]). In absence of any orientational order, the shape of the 
interface is given by the classic catenoid, parametrized as: 



X = a cosh ^— ^ cos , ^ = a cosh ^— ^ sin ^ , 



z = u 



where < (/> < 27r is the usual polar angle and —h< 
u <h. The mean curvature is everywhere zero while the 



Gaussian curvature is given by i^^ = — sech"^ (ix/a) /o?. 
The length a is obtained by setting the polar radius to 
R when z = ±h. This yields R = acosh{h/a). It is well 
known that the solution of this transcendental equation 
exists only if h < R/ mina;{cosh(x)/x} :^ 0.663 i?, while 
for h/R> 0.663 the catenoid does not exist. In the {u^ (j)) 
coordinates the spin connection has components: = 
and A(f) = — tanh('u/a). 

For large orientational stiffness we expect the inter- 
face bounded by the two circular rims to be a surface of 
revolution of constant negative Gaussian curvature. A 
surface with this properties can be identified in the pseu- 
dospherical hyperboloid whose parametric form is 
given by: 
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with < < 27r and -/7 < < /7. Here E{(l)\m) = 
dt{l — msin^t)^/^, with \m\ < 1, is the incomplete 
elliptic integral of second kind. The coordinate u is lim- 
ited in magnitude \u\ < aasinh(a/6) and the surface has 
constant negative Gaussian curvature K = —1/a^. This 
immediately fixes the value of the constant a = (/^/<j)^/^. 
The length 6, on the other hand, is found by requesting 
z = when u = ±U = ±a acosh(i?/6). This gives: 
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(8) 



The solution of this equation must satisfy a > R > b 
and R^ < + b^ as the consequence of the bounds 
on the coordinate u and the elliptic modulus. The 
components of the spin connection are: Au = and 
^0 = sinh(ii/a). 

Now, because both the catenoid and the pseudospher- 
ical hyperboloid are azimuthally symmetric, all the rel- 
evant quantities will be independent on (/), in particu- 
lar = 0{u) and cp = (p{u). This latter property im- 
plies that 6 is constant throughout the interface since: 
duO = Au -\- ^t^(f)^ = 0. Thus the orientational energy 
density simplifies to |Ap = g^^A^ {gij being the metric 
tensor). The total energies can be readily calculated in 
the form: 
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— (X - atan x) H 

a K, 



, (9a) 



(9b) 



where x = V {R/bY — 1. As a consistency check, one 
could also verify that the geometric potential ip is given 
for both surfaces by the simple function: ip = log(r/i?), 
with r the usual polar distance in the xt/— plane. Sub- 
stituting in Eq. ([5| and integrating, promptly gives Eqs. 
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FIG. 2: Phase-diagram obtained by comparing the energies in 
Eqs. Q. For gF? j k > 1 the equiUbrium shape is a catenoid. 
aR^/hv < 1, on the other hand, the energy balance is domi- 
nated by the orientational stiffness of the ordered phase and 
the interface is a surface of constant negative Gaussian cur- 
vature K — —G I K. 

(|9|. A comparison of the energies ([9| leads to the dia- 
gram of Fig. [2] in the plane [gR? j K^hj K). As expected 
for (jR^ /n > 1, surface tension dominates and the inter- 
face has the classic shape of a catenoid. For cjR^ /n < 1, 
on the other hand, the energy balance is dominated by 
the orientational stiffness of the ordered phase and the 
interface adapts to the in-plane orientational order by 
forming a surface of constant negative curvature. The 
stability of this purely hyperbolic interface is enhanced 
when the aspect ratio h/R is closer to the upper bound 
h/R - 0.663. 

A check of the orders magnitude is inevitable at this 
point. Standard liquid crystals have ~ /c^T ~ 10~^^ 
erg at room temperature, whereas the typical surface ten- 
sion of liquid crystal interface is of order a ~ 10 erg/cm^. 
Thus, having aR^ /n ~ 1 would require a rim of radius 
R < 1 nm, which is smaller than a typical film thick- 
ness. Block copolymer membranes, as those considered in 
Ref. [5], are more suitable candidates to observe the hy- 
perbolic morphology described here, as the orientational 
stiffness is of order 60 /c^T. In actomyosin compounds, 
the contractile forces generated by the molecular motors 
have been recently demonstrated to give rise to an effec- 
tive surface tension p!5]. Due to their elongated shape 
as well as the cross-linking operated by the myosin mo- 
tors and, possibly, by other actin-binding-proteins, actin 
filaments can form a nematic phase at room tempera- 
ture. Since the forces exerted by the motors, of order of 
pico-Newtons, are in this case responsible for both the 
surface tension and the alignment, one can expect their 
effect to be comparable in magnitude so that the balance 
(jR^ / n ~ 1 could be achieved for reasonable system sizes. 

Hexatic suspensions obtained by trapping colloidal 



particles at the interface between two liquids, represent 
an especially promising candidate for an experimental re- 
alization the phenomenon discussed in this Letter. In 
this case the orientational stiffness can be finely con- 
trolled by moving across the liquid/hexatic/solid phase 
diagram. As predicted in the theory of two-dimensional 
melting [16 the orientational stiffness has a universal 
value n/ksT = 72 /n at the liquid/ hexatic transition 
and diverges exponentially at the hexatic/solid transi- 
tion. This was brilliantly verified by Keim et al. p/7j in 
experiments with superparamagnetic colloidal particles 
with tunable interaction strength at the water/air inter- 
face. Thus, by approaching the hexatic/solid transition, 
one could in principle make the orientational stiffness ar- 
bitrarily large (~ 10^ ksT in practice). In turn, the 
interfacial tension between two immiscible liquid phases 
can be made arbitrarily small as the critical point for 
phase separation is approached. 

In conclusion, I showed that surfaces of constant nega- 
tive Gaussian curvature can be realized in fluid interfaces 
equipped by an orientational ordered phase as a conse- 
quence of the interplay between surface tension and ori- 
entational elasticity. A simple example of this behavior 
was presented in the case of defect-free tubular interfaces. 
An experimental realization of this phenomenon could be 
achieved in hexatic colloidal suspensions and remains a 
challenge for the future. 
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